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Abstract 



A theory of pump-probe spectroscopy is developed in which optical fields 
drive two-photon Raman transitions between ground states of an ensemble of 
three-level atoms. Effects related to the recoil the atoms undergo as a result 
of their interactions with the fields are fully accounted for in this theory. 
The linear absorption coefficient of a weak probe field in the presence of two 
pump fields of arbitrary strength is calculated. For subrecoil cooled atoms, 
the spectrum consists of eight absorption lines and eight emission lines. In the 
limit that xi <C X2> where xi an d Xi are the Rabi frequencies of the two pump 
fields, one recovers the absorption spectrum for a probe field interacting with 
an effective two-level atom in the presence of a single pump field. However 
when xi ^ X2, new interference effects arise that allow one to selectively turn 
on and off some of these recoil induced resonances. 



I. INTRODUCTION 

Recent advances in laser cooling, atom optics, and Bose-Einstein condensation have 
underlined the role played by atomic recoil in atom-field interactions. A measure of the 
importance of recoil effects is the recoil frequency, = hk 2 /2M, associated with the 
absorption, emission or scattering of radiation of wavelength A = 2ir/k by an atom of 
mass M. Once this quantity becomes greater than or comparable to decay rates or Doppler 
widths that characterize the spectral response of atoms, recoil can lead to new features in 
absorption or emission line shapes. One class of such phenomena has been termed recoil- 
induced resonances (RIR) [1-8], which occur when a weak probe and strong pump field 
simultaneously drive a given atomic transition. Interesting in their own right, the RIR have 
been used to determine the velocity distribution of laser-cooled atoms 0, as a probe of 
Bose-Einstein condensates || |10[| , and in a feedback mechanism in stochastic cooling flTT 



Related to the RIR is the so-called collective atomic recoil laser (CARL), which operates on 



similar principles but in a somewhat different parameter range |T2j . Both the RIR and CARL 



represent new diagnostic probes of cold-atom sysytems. Recently, Moore and Meystre |I3 
proposed that CARL be used to entangle optical and matter fields, as well as to entangle 
different modes of the condensate excited by the optical fields. Our discussion is limited 
to situations in which the collective effects associated with CARL can be neglected. In 
this paper, we combine RIR with ground state spectroscopy ]14|] to obtain qualitatively new 
features in the probe absorption spectrum. 



The scheme we adopt is based on the model developed in [|14[| involving Raman transi- 
tions, but for which all effects associated with atomic recoil were ignored. In that work a 
new type of interference was discovered, allowing one to selectively turn on and off certain 
lines in the absorption-emission spectrum by controlling the ratio of the Rabi frequencies of 
the two fields that comprise the two-photon pump field. Interference in a dressed state basis 
occurs between pathways involving the probe field and each of the two pump fields. Since 
the two pump fields impart different recoil momenta to the atoms, it is not at all obvious 
that the interference persists when recoil splittings are resolved. Part of the motivation 
for our calculation is to examine this question. In addition, we show that the interference 
persists even if the pump fields are in quantized, Fock states. 

The probe absorption spectrum consists of as many as eight absorption-emission doublets 
which are fully resolvable if ujk > 7, where 7 is some effective ground state lifetime. This 
is in contrast to the RIR spectrum on dipole allowed optical transitions |l|], where at most 
one absorption emission doublet is resolvable if uo^ < 7 e , where 7 e is an excited state decay 
rate. One might question the need to increase the number of recoil doublets in the probe 
spectrum, since a single doublet can be used to probe recoil effects. We show that the 
additional recoil structure reflects interesting quantum dynamics of the combined atom-field 
system, as well as providing some new applications. 

In Sect. II, a model is developed for the interaction of the atoms with the pump fields; 
dressed states of the atom plus pump fields are defined. In Sect. Ill, the interaction with the 
probe field is introduced and the dressed state picture is used to obtain the probe absorption 



spectrum in the secular limit. In Sect. IV we discuss the results and possible applications. 
Nonsecular contributions to the absorption coefficient are calculated in an Appendix. 



II. EQUATIONS OF MOTION 

The atom-field system is indicated schematically in Fig.l. Ground state levels |1 > and 
1 2 > are pumped incoherently with rates Ai(p) and A2(p), respectively, and the population 
of both states decay with rate 7. If states |1 > and |2 > represent stable ground states of 
the atom, then the pumping rates and decay rate constitute a simple model for atoms that 
enter and leave the interaction volume. The ground to excited state transition frequencies 
are denoted by u e j (j = 1,2). The pump fields 1 and 2, which constitute the two-photon 
pump field, are denoted by the coupling constants g\ and g 2 , respectively. Pump field 1 
couples only state |1 > and excited state |e >, while pump 2 couples only states |2 > 
and |e >. The pump fields have frequencies Qi and f2 2 , and propagation vectors kx and 
k 2 , respectively. In this section, equations are derived for the atom-pump field interaction, 
neglecting the incoherent pumping and decay of the ground state levels. In the following 
section, the interaction of the atoms with the probe field, E p , which couples only states |1 > 
and |e >, is incorporated into the calculation, as are the incoherent pumping rates Ai(p) 
and A 2 (p), and the ground state decay rate 7. 



In contrast to |L4j but as in 0, we use a quantized description of the pump fields. If the 



pump field detunings, 

A 2 = fi 2 - oo e2 ps Ai = O x - u el = A, (2) 

are sufficiently large such that |xi,2/A| 2 <C 1 and |7 e /A| <C 1, where the xi,2 are defined 
below and 7 e is the excited state decay rate, it is possible to adiabatically eliminate the 
excited state to arrive at an effective Hamiltonian involving only states |1 > and |2 > which 
is of the form 
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H = H a + H r + H ar ; 



2 2 

(huji + — )|l,p >< l,p| + {hu 2 + — )|2,p >< 2,p| 

2M 2M 



ft — \2, p >< 1, p|a 2 ai + n — 11, p >< 2, p\a[a 2 



(3a) 
(3b) 

(3c) 
(3d) 



H r = hVL\a\ai + M^l^i 

H ar = / ] 

p 

where H a is the Hamiltonian for the atom in which the center-of-mass momentum p has 
been quantized using periodic boundary conditions in a volume V (assuming that the atoms 
are free and not subject to some trapping potential), H r is the free field Hamiltonian for the 
two pump fields, H ar represents the interaction of the two-photon pump field with an atom, 
in the rotating- wave approximation, and 



Bare states, \j, p; n 1; n 2 >, are defined to be eigenstates of H a + H r , where j — 1,2 labels the 
internal state of the atom and ri\ and n 2 are the number of photons in pump fields 1 and 2, 
respectively. A term in the Hamiltonian corresponding to the light shifts of the ground state 



levels, Hi s = 2~2p ^~5~|1?P >< l>p| a I a i + ^^^|2,P >< 2,-p\a\a 2 , has been omitted in 
Eq. (|3aD ; such light shifts can be included by a redefinition of the ground state frequencies, 
^2 + (a 2 a 2j —> ^2 and u 1 + (a^aA — ► uj x . 

The matrix elements of the operator e* k R in the momentum-state basis are 



92 1 



< p|e ik ' R |p' >=< p|p' + hk >= 5 P:p/+hk . 



This allows one to rewrite the interaction term as 



h 



9i9 2 



2, p + hk 12 >< 1, p|<4ai + h-^-^\l, p >< 2, p + hk 12 \a\a 2 



(4) 



(6) 



The Hamiltonian H results in an infinite ladder of decoupled two state manifolds (p, rii,n 2 ) 
involving the states |l,p;ni,7i2 > and |2,p + ftk 12 ;ni — l,n 2 + 1 >. The Hamiltonian for 
the manifold (p, ni,n 2 ) is 



if (p, m, n 2 ) = e(p, m, n 2 )I + h 



f -~S(p)/2 G* N 
v G *(p)/2 ; 



(7) 



where I is the identity matrix 



G 



X1X2 
A 



\G\e 



h 

e(p,ni,n 2 ) = fr(ni^i + n 2 Q 2 ) + ^(^p + ^p+^ 12 + ^1 + ^2 + ^2 - fii] 



<Kp) = ^12 - VKki2 - 



pk 



12 



M 



5 U = Ai - A 2 = - fi 2 ) - ^21, 



^ P = ^7, and \2 = W™2 + 1, Xi = 9\^\- 

The dressed states are defined to be the eigenstates of the matrix in Eq.(^), with energies 

771 / \ / x , ^ab(p) 
E b ,a\P) = £{P,ni,n 2 ) ± ; 



and associated eigenkets |L5 
/ 



where 



and 



V 



\A > 
\B > 



J 



T*(p) 



^ |2, p + /lki 2 ; ni — 1, n 2 + 1 > ^ 



T(p) 



oi<t>d/2 



V 



cos [0(p)] 
sin [0(p)] 



|l,p;ni,n 2 > 



_ p -*0d/2 



sin[0(p)] 



COS [0(p) 



COS [0(p)] 



*(p) \ 



1/2 



(9) 



(10) 



(12) 



The value of 0(p) is restricted to the range < 0(p) < 7r/4 for <5(p) > and 7r/4 < 6(p) < 
7r/2for 5(p) < 0. For 0(p) ~ (5(p) > and 5(p)/|G| > 1), \A >~ |2, p+£k 12 ; m-1, n 2 + 
1 > while for 0(p) ~ vr/2 (5(p) < and <5(p) > 1), |-B Q >~ |2, p+ftk 12 ; m-1, n 2 + l >. 
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III. PROBE FIELD ABSORPTION IN THE SECULAR LIMIT 



The effect of the probe field is to induce transitions between states in different manifolds. 
As is customary in dressed atom approaches, the probe is treated as a classical field, 

E(R,t) = heE^- 11 -***) + c.c, (13) 

where e is a unit polarization vector. For our problem, however, this choice represents 
a hybrid approach, since two-quantum processes involving the probe field and either of the 
pump fields mix classical and quantized fields. Although the probe field is treated classically, 
its effect on the momentum of the states must be accounted for explicitly. If the probe field 
detuning on the 1 — > e transition is sufficiently large to be consistent with the adiabatic 
elimination of the excited state, all transitions involving the probe occur via two-quantum 
transitions involving the probe field and either of the pump fields. In the bare state basis, 
starting from the (p, 711,712) manifold, probe field absorption corresponds to transitions 
|l,p;7ii,7i 2 >^ 1 1, p + ftkpijTii + l,n 2 > or |l,p;7ii,7i 2 >^ |2, p + hk p2 ; Tii, n 2 + 1 >, 
where the second photon is emitted into either the pump 1 or pump 2 modes, respectively. 
Similarly, probe gain corresponds to transitions |1, p; ni, n 2 >^ 1 1, p — Kk p i, 71 1 — 1, n 2 > or 
|2, p + Kki 2 ; rii — 1, n 2 + 1 >^ |1, p — Kk p i, n x — l,n 2 >. This is illustrated in Fig. 2(a). 

This picture of probe field absorption or emission allows us to reach an important con- 
clusion concerning interference between pathways involving both pump fields. The two 
absorption processes shown in Fig. 2(a) involve different final states and do not interfere, 
reinforcing the possibility mentioned in the Introduction that interference may be suppressed 
when recoil is taken into account. However, both final states belong to the same manifold - 
the (p + Kk p i, rii + 1, n 2 ) manifold. As such, when these states are dressed by the two-photon 
pump field, each state in the final state manifold will be coupled to the initial state by two 
separate pathways involving the probe field and each of the pump fields. This implies that 
terms in the probe absorption depending on the simultaneous presence of both pump fields 
will exhibit interference effects. In the case of probe gain, it is the initial states that differ, 
but the overall conclusion remains unchanged. 
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The dressed state approach provides a convenient and relatively easy method for obtain- 
ing the probe absorption spectrum in the secular limit, where the frequency separation of 
the dressed states in a given doublet is much larger than the ground state decay rate 7. [It 
is assumed from this point onward that 7 e (xi. 2 /A) 2 <C 7, implying that the dressed states 
decay with rate 7 One need calculate only the dressed state energies and the transition 
matrix elements to obtain the spectrum. A more detailed treatment of the problem, allowing 
one to calculate non-secular contributions, is presented in Appendix B. 

It is straightforward to generalize the dressed states defined in the previous section 
to include the two manifolds coupled to the initial manifold by the probe. The 0, 
1, and 2 manifolds refer to (p, n 1 ,n 2 ) = {|1, p; ni, n 2 >, |2, p + ftki 2 ; rii — 1, n 2 + 1 >}, 
(p + frkpi,?2i + l,n 2 ) = {|1, p + Kk pl ; rii + 1,^2 >, |2, p + h(k 12 + k pi ); n x ,n 2 + 1 >} , and 
(p - hkp^ni - l,n 2 ) = {|l,p - hkp^ni - l,n 2 >, |2,p + h(k 12 - kpi);ni - 2,n 2 + 1 >}, 
respectively. Taking the central energy of the initial manifold (p, rii,n 2 ) arbitrarily equal to 
zero, one finds that the dressed state energies are given by 

E% = ^ 10 (p) ± ^S(p); (14) 
e a]b = ^20 (p) ± ]^wf B {p), 



where 



^o(p) = <Hp) = £12 - ^hk 



ou%(p) = ^A\G\ 2 + ^V)\ (15) 
P ■ ki2 



M 



K ( \ It 1 fcl \ x (p + fikpi) ■ ki2 . , 

di{p) = d{p + hk pl ) = b 12 - u h]il2 — ; (16) 



8 2 (p) = S(p - hkpt) = 5 12 - u hkl 



M 

[p - fikpi) ■ kj 2 
M 



^io(p) = e(p + fikpi, n x + 1, n 2 ) - e(p, m, n 2 ) = h\Vt l + u hkpl + kp j^ P + —^jf^ J 5 ( 17a ) 
^2o(p) = e(p - Kkpt, ni - 1, n 2 ) - e(p, m, n 2 ) = ( -fii + ^ kpl - ^ P - ^\ kl 2 ) , (17b) 



and it has also been assumed that n\,n2 ^> 1 such that G^ -* = g\gi-\ln\ (n 2 + 1)/A m 
= gtg*y/{jix + 1) (n 2 + 1)/A « G( 2 ) = ^ 2V /(m - 1) (n 2 + 1)/A = G. The dressed 
state angles are given by 



COS [0j(p)] 



and dressed state kets are defined by 



1 + 



ft(p) 
^ab(p) 



1/2 



(19) 



/ 

V 



Si > 



\ 



/ 



T*(p) 



2, p + h [k 12 -(-l) 1 (1 - 5 ifi ) kpi] ;m - 1 - (-1)* (1 - , n 2 + 1 > 
|1, p— (— 1)* (1 -5i,o)ftkpi; w i - (! - fc,o) ,w 2 > 



\ 



(24) 



where 



T*(p) 



e i0 d /2 cos [Q.fa)] - e -i*d/2 S i n [#.( p )] 

e i0 d /2 gin . ( p )] e -i<t> d /2 cog ^. (- p )] 



(22) 



The absorption coefficient, a, is proportional to the rate at which energy is absorbed 
from (a > 0) or emitted into (a < 0) the probe field. Absorption corresponds to transitions 
from initial dressed states / = A ,B to final state dressed states J = Ai,Bx, while emis- 
sion corresponds to transitions from initial dressed states I = A , B to final state dressed 
states J = A 2l B 2 - For a given transition, the contribution to the absorption coefficient is 
proportional to NVhfl p '-ypjj/\E p \ 2 , where N is the atom density and pjj is the steady state 
population in state J owing to the / — > J transition. The final state population pjj(p) is 
equal to [A/(p)/7] | < J|V^|7 > \ 2 L U , where A/(p) is the pumping rate for initial dressed 
state / [16|], < J|lp|/ > is a matrix element for the I — J transition and Ljj is a Lorentzian 
having width 7, centered at the I- J transition frequency. The transition frequencies A/j(p) 



may be directly read from Fig. 2(b) or obtained from Eq. (13). Consequently, it is neces- 
sary to specify A/(p) and to calculate the transition matrix elements < J|V^|J > in order 
to arrive at an expression for the probe absorption. 

In the bare state representation, the pumping matrix for the intitial state manifold is 
taken to be of the form 



9 



A(P) = 



\ 



pump 



PnlP) o 
^ p 22 (p+ftki 2 ) J 

which, when converted to the dressed state basis becomes 
/ 

Mp) = 



Ai(p) 
A 2 (p + hk 12 ) 



\ 



A A (p) Aab(p) 
Aab(p) A B (p) 

A 2 (p + Hk 12 ) cos 2 [0 o (p)] + Ai(p) sin 2 [0 o (p)] \ [A 2 (p + hk 12 ) - A^p)] sin [20 o (p)] 

| [A 2 (p + fik 12 ) - Ai(p)] sin [20 o (p)] Ai(p) cos 2 [0 o (p)] + A 2 (p + fik 12 ) sin 2 [0„(p)] 



In the secular limit, the off-diagonal pumping terms can be neglected since they give rise to 
terms of order 7/^5 (p) <C 1. Moreover, for the present we will set A 2 (p + Kk 12 ) =0 and 
generalize the results to nonvanishing A 2 (p + Kk 12 ) in the next section. Thus, we take a 
pumping matrix of the form 

/ 

Ad (P) = 



A^sin 2 [0 O (p) 







V 



Ax (p) cos 2 [9 (p)] 

We now turn our attention to the transition matrix element < J|lp|J > . As an ex- 



ample consider < A^V^^o >, which is the amplitude for the transition \A \A 1 > 
involving probe absorption and pump 1 or pump 2 emission. This transition is illustrated 
in Fig. 2(c). The probe couples only to the |l,p;ni,n 2 > part of \A Q >, leading to a factor 
_ e *0d/2 g j n [6» (p)]. The absorption of the probe is followed by emission into pump 2 taking 
the atom to the |2, p + h(k±2 + k p i); n\, n 2 + 1 > component of \A 1 > and leading to a factor 
e -«<?V2 cos [6' 1 (p)] or emission into pump 1 taking the atom to the the |1, p + ftk pl ; ni + 1, n 2 > 
component and leading to a factor — e i<?id / 2 sin [0i(p)]. The coupling strengths for these two, 
two-photon transitions are G\ and G 2 , respectively, where 



(25) 



The two processes add coherently and one finds 

< A 1 \V\A >=h{G* 2 e^co S [^(p)] - G>n [^(p)]} {- sin [0 o (p)]} e "«H 



10 



Other matrix elements are calculated in a similar manner. For probe gain, pump fields 1 and 
2 couple to the |1, p; rii, n 2 > and |2, p+Kk 12 ; n 1 — 1, n 2 + l > components of the dressed states 
in the manifold, respectively, while the probe field couples to the |1, p — hk p i, rii — 1, n 2 > 
component of the dressed states in the 2 manifold. Explicit expressions for the matrix 
elements are given in Appendix A. 

Combining all transitions and summing over p, one finds an absorption coefficient pro- 
portional to 

7 2 A/(p) ( | < J\V P \I > | 2 v- \<J\V P \I>\ 2 

(2 

For a sub-recoiled cooled vapor, we can set A^p) = Ai5 p , , such that A j4 (p) = A A 5 Pt0 and 
A s (p) = A B 5 pfi , where 

A^ = Aisin 2 ^ ; Ab = Aicos 2 ^, 

^ = 0j(O) [and, for future reference, A/j = A/j(0), u;^ = a;^(0), etc.] and 5 p0 is a 
kronecker delta. In this limit, one finds the absorption coefficient in the secular limit to be 
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a 



IAI 



^ sin ± _ I cos fl^ ^ sin 2 ^ LaoAi ( A /) + ^ cos 2 ^ LsoAi (A') 

Immsfl, 4- -sinfl^ f — sin 2 Q n L ... „ .(A') 4 — cos 2 n /", D . „. CA'^ 



(A') + ^cos 2 ^ L BoBl (A') 
1 / V7 7 

- ^77 sin 6> - i cos 6» ^ ^ (sin 2 6> 2 La a 2 (A') + cos 2 9 2 L AoB . 2 (A')) 

1 . \ Ar / r> „ _ , . ,, r, „ _ , . ,, N 

V 



4 ( ^77 cos 61 4 - sin ^1 J ( — sin 2 L AoBl 
\ V J V7 

/ 1 \ 2 A . 



(1 \ 2 A 

V>?7 cos 6 + - sin # J (sin 2 2 £b o a 2 (A') + cos 2 6 2 L Bo b 2 (A')) 



7 2 


7 2 


2 

+ f 


(A'-a; nkpl + I(4°) 
7 2 


2 

+ f 


[A> - u mpl - \{uf B + 
7 2 


2 

+ f 


(A'-a; nkpl + i(4°i-a;«)) 
7 2 


2 

+ f 


(A' + u mpl + \{uf B - uf B )) 
7 2 


2 

+ f 


(A' 4-^ + 1(4°) +uf B )) 
7 2 


2 

+ 7 2 


(A' + u hKl -\{uf B + ujf B )) 
7 2 


2 

+ 7 2 


(A' + co hkpl - \{uf B - u% )) 


2 

+ f 



LbqA\ (A') 
Lb q Bi (A') 

^ A 2 (A') 
^ B 2 (A') 
^oA 2 (A') 
^B B 2 (A') 

where 

^ = A/|A| 

is the sign of the detuning for each of the three fields, 



I (27a) 
(27b) 

(27c) 

(27d) 

(27e) 

(27f) 

(27g) 

(27h) 

(27i) 
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, hk pl ■ kia 

A — dpi — —Mr- , 

8 p i = Q p — Oi ; (5 p 2 — — ^2 — ^2i) 



k p Nd( 

a = 



2 



2he 'j ' 

and di e is a bare state dipole moment matrix element for the |1) — > |e) transition. 

The first four resonances in ( — J correspond to probe absorption while the last four 

V a ° / sec 

correspond to probe gain. The spectrum is shown in Fig. 3. The line widths of all the 
resonances equal 7; consequently, the recoil induced resonances in ground state spectroscopy 
are fully resolved if Uk > 7- The secular contribution does not vanish in the limit that 
\5(p — 0)/G\ < 1, even though = A B in this limit. As long as the recoil frequency is 
larger than 7 and the Doppler width associated with the two-photon pump transition, the 
absorption and emission contributions to the probe response do not cancel one another. 

The most significant feature of (^^j is that the line strengths involve factors such as 
[ipr] sin d\ — ^ cos which allows one to manipulate the strength of the lines by controlling 
the sign of the field detuning and the ratio of the pump field amplitudes. These factors are an 
indication of interference between the two "two-photon probe" fields which can both lead to 
absorption or gain in E p . Because the two lines in the doublets have different strengths, one 
can adjust ip and 77 to turn off one of the lines. For example, the absorption doublet, La a 1 
and LbqBx, consists of the lines at A' = ^k pl + 1(^ab~ u ab) an d - |(^ab ^^ab) witn 
strengths ~ (ipr) sin 9i — ^ cos O^j and ~ (^prj cos 61 + i sin , respectively. Choosing 
if; — +1 and rf = cot9 1 turns "off" the first line while if; — — 1 and rf = tanOi turns "off" 
the second line. This is shown in Fig. 4. When \G\ and 812 are much larger than any of 
the recoil terms, #0 ~ $1 ~ $2 and the emission lines are also turned "off". Consequently, 
by choosing if; — +1 and rf = cot#i to turn off the La {) a x and Lb a 1 absorption lines, the 
La q a 2 and La b 2 emission lines are also turned off. 

A particularly interesting case occurs when k x k 2 so that uu^b = u ab = u ab an d $0 — 
61 = 9 2 . This would correspond to a two-photon pump field which imparts no momentum 
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to the atoms so that the recoil splitting in the absorption spectrum can be attributed solely 
to the recoil due to the probe field acting with either of the pump fields, k p2 = k pl . In 
this case, the line La a 1 (A') is degenerate with Lb^b^ (A') and La a 2 (A') is degenerate 
with Lb b 2 (A') . Consequently, the spectrum consists of three absorption-emission doublets 
centered at A' = 0, +oo^ B , —oj^b- Moreover, the lines within each doublet are split by 2ujf^ 
which is independent of the strength or detuning of the pump fields. 

When the effects of atomic recoil are neglected by setting all recoil momenta to zero in 
(j^J > one obtains the same absorption spectrum given in [14|. In the limit that G\ = 0, one 
recovers a simple, recoil shifted Raman spectrum. In the limit that G 2 = 0, one recovers 
the central, secular components of the pump-probe spectrum associated with a single, two- 
level optical transition jj. In the limit that r\ <C 1, while G/7 ^> 1 remains constant, the 
absorption spectrum mirrors that for the pump-probe spectrum associated with a two-level 
optical transition [|J. 



IV. DISCUSSION 

For a subrecoil cooled vapor, pumping to state 2 at a rate A 2 (p + Kk 12 ) = A 2 5 p+ ^k 12! o 
doubles the number of absorption and emission lines in the probe spectrum, but does not 
result in any qualitatively new features. There will be an additional contribution to Eq. 
(H) in which Q { = 6^(0) is replaced by 0* = 0j(-ftk 12 ), A u by A u = A u (-hk 12 ), u % AB 
by &ab = ^abC - ^12), by A 2 cos 2 ^o and A# by A 2 sin 2 ^o- The absorption coefficient 
contains sixteen lines in all. The eight new spectral components display the same properties 
as the original eight but are displaced by an amount ~ ^pl, F i gure 5 shows the secular 
absorption spectrum with all sixteen components when A 2 = Ai. In the absence of recoil, 
the secular absorption coefficient vanishes when A 2 = Ai ||14|. However, when recoil is 
included, the absorption and emission contributions to the probe response do not cancel one 
another when A 2 = A 1; provided the recoil frequency is larger than 7 and the Doppler width 
associated with the two-photon pump transition. 
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In Ref. |14] interference effects similar to those in Eq. ( |27a|) were found. We have shown 



that the interference is independent of the field statistics and persisits even when recoil 
induced resonances are resolved. 

The RIR offer several possibilities for applications. The existence of a number of tunable, 
well-resolved, gain peaks allows one to envision experiments in which lasing occurs at one 
or more of the gain positions. These peaks could be adjusted to coincide with modes of a 
ring cavity, for example. By sweeping the ratio of the two pump field intensities, one has 
a mechanism for modifying the probe gain or absorption, to the point of total suppression. 
The central frequency of the absorption-emission doublets can be controlled via pump field 
strength and detuning. For parallel pump fields, the frequency separation of the absorption 
emission doublets can be as large as 80;^ and is independent of pump field strength and 
detuning. With line widths approaching 1 Hz or less, the probe spectrum can be used to 
measure the recoil frequency to a precision of order 10~ 7 ; this precision can be increased if 
the two-photon pump field is replaced by a pair of counterpropagating, two-photon pump 
fields []17| . The narrow resonances can also be used in schemes for obtaining "slow light" 

Pump-probe spectroscopy of Bose-Einstein condensates represents an interesting appli- 
cation of the ideas presented here. Bragg spectroscopy has recently been demonstrated in 
condensates |J [[K]] as well as the stimulated generation of matter waves in a condensate 
by Rayleigh scattering |19|| . Pump- probe spectroscopy using electronic excited states would 



be unfeasible in condensates since there are no stable trapped condensates with electronic 
states which may be populated. Currently, the only multi-component condensates consist of 
two hyperfine states in 87 Kb J21] and the Zeeman states in the F = 1 manifold of optically 
trapped Na [|22[ . Consequently, pump-probe spectroscopy would necessarily involve Raman 
transitions between stable ground states in the manner proposed here. The RIR spectrum 
of a weakly interacting Bose condensate should yield information about the spectrum of 
elementary excitations in the condensate which for small momenta have a linear dispersion 
relation while, for large momenta, have a quadratic dispersion similar to that of free atoms, 
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but with a shift due to the mean-field interactions in the condensate. In addition, the line- 
widths of the RIR spectrum should be given by the zero point motion of the condensate in 
the trapping potential provided (^^) k p i >2 > 7 and Aa; is the size of the condensate. How- 
ever, the direct application of the results presented here to a condensate would be erroneous 
since a correct calculation of the RIR spectrum would have to account for the mean-field 
interactions between the atoms. This will be pursued in future work. 
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VI. APPENDIX A - MATRIX ELEMENTS OF \ ID 

The interaction with the probe field may be expressed in terms of an effective two-photon 
interaction Hamiltonian similar to that in Eq. (|3d|), 

Vi = hJ2 O^e-** |2, p + hk p2 ) (1, P | 4 + Ml e -«v |i ; p + hkpl) (i, p| fl t + h A 

p 

(28) 

where \p — ^kdeiEp is the Rabi frequency for the probe field and d e \ = < e|d-e]l > is a 
dipole matrix element. The matrix representation of V/ with respect to the bare state basis 
has the following nonvanishing elements 

(1, p; m, n 2 \ V/ |2, p + Kk 12 + Kk pl ; n 1 ,n 2 + l) = hG 2 e iQpt ; 

(l,p;n 1 ,n 2 \V I \l,p + hk pl ;n 1 + l,n 2 ) = hG 1 e i ^ t ] (29) 
(1, p-ftkpi; m - 1, n 2 \ V/ |2, p + Kk 12 ; n x - 1, n 2 + 1) = hG 2 e iQpt ; 

(l,p-fikpi;m - l,n 2 \ Vj |l,p;ni,n 2 ) = hde^ 1 ; 
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and the hermitian conjugates of Eqs. (29). In Eqs. (29), G\ and G 2 are two-photon probe 
Rabi frequencies defined as G 1 = x* v X\l 'A; G 2 = X* p X2j 'A. Note that the coupling of the 1 
and 2 manifolds to manifolds other than the manifold has been ignored. 

The matrix V Id represents the interaction with the probe field in the dressed state basis 
and is defined as 

V Id = TV/T 1 "; (30) 
where T is given by the block diagonal matrix 

T = ^[T 1 (p),T (p),T 2 (p)]. (31) 
The Rabi frequencies Xpi Xi> an d X2 ma Y be expressed as 

Xv = \Xp\e i(p ; Xi = IXile^ 1 ; X2 = \X2\e^; (32) 
so that <fi d = 4>2 — 0i + f (1 — VO since G = \G\e l<t>d . The matrix elements are: 



<A l 


v, 


A 


>= 




h) e - 


-^'sin0 o ( 


— \G 2 \ cos^! + ip\G\\ sin^i); 


(33a) 


<A, 


v, 


B 


>= 




h) e - 


- <n **cos0 o 


' \G 2 \ cos 61 — i)\Gi \ sin#i); 


(33b) 


<B 1 


v 7 


A 


>= 




h) e - 


sin O ( 


— i^\G]\ cos^i — \G 2 \ sin^x); 


(33c) 


<B 1 




Bo 


>= 




h) e - 


cos e 


[tp \Gi\ cos^i + \G 2 \ sin^); 


(33d) 


< A 


v, 


A 2 


>= 


he^-' 


h) e - 


-^*sin 9 2 ( 


— \G 2 \ cos # + -0 Gi sin 9 ); 


(33e) 


< A 


v, 


B 2 


>= 




h) e - 


- <n "*COS0 2 


[\G 2 \ cos6 - 4>\Gi\ sin6* ); 


(33f) 


< B 


v 7 


A 2 


>= 




h) e - 


- <n **sin0 2 ( 


— if)\Gi\ cos 6*0 — G2 sin^o); 


(33g) 


< B 




B 2 


>= 




h) e - 


cos e 2 


[ip \Gi \ cos# + G2 sin# )- 


(33h) 



The other elements follow from the hermiticity of \i d . 

VII. APPENDIX B 

In this Appendix the absorption coefficient for the probe field is calculated, without 
making the secular approximation. The absorption coefficient, a, and index change, An, 
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arise from the imaginary and real parts of the macroscopic polarization in the Maxwell-Bloch 
equations for the probe field. They are given by the expressions 



v 

At, 



a = a Im ( ^-V ) ; (34a) 



An = -aok- 1 Re [j^ V J 5 ( 34b ) 
kpNdf 

ao = ~2H^- (34C) 
where V is the volume and p' le is the part of the bare state density matrix element pi e (R, t) 
which is proportional to e -*( k p- R -^p*) ; which we denote by p' le (R, t). 

Before proceeding, we note that in this Appendix all summations over momentum states 
have been converted to an integration over a continuum of states via the standard substi- 
tution J2 P ~^ (2nh) 3 I ^ 3 ^' ^ ne coe ffi c i en t Pie i s related to the coherence in position space, 
p'i e (R, t), and the momentum space density matrix elements, pi e (p, p'; t) = p'i e (p, p'; t)e tnpt , 
by 

p'ie(R^) = p' le e-^- R -^ (35) 
= j^~ 3 J J d 3 pd s p'p' le (p, p'; t)e^e^~^ h 5(p - p' + ftk p ). (36) 

The coherence, p' le (p, p'',t)e pt , has been written in the Schrodinger representation and is 
obtained from the density matrix for the atom plus pump fields by tracing over the number 
of photons in the pump fields, 



P 



' le (p, p'; t) = e-^'e— ' ]T p( e (p, m, n 2 ; p', m, n 2 ; t). (37) 



nx,n2 



where p{ e (p, rii, ni\ p', n^, n^; t) is in the interaction representation with respect to the inter- 
nal energy levels and pump fields. One cannot derive a differential equation for p le (p, p'',t) 
starting from the original Hamiltonian in Eq. ([3a]) since the excited state has been adiabat- 
ically eliminated from the effective Hamiltonian. However, by reintroducing an interaction 
term 



H af = h 



0i |e> (1| e iklR + g 2 \e) (2| e ik2 ' R + Xp \e) (1| e< kl ' R - n * f ) 



+ h.c. 



into the Hamiltonian and writing 



p' le (p,p';£) = ^ L~^ e -^ PufarHMlf, 

V ni,n 2 

one finds the equation of motion for p{ e (p, nx^n^ p', n[, n 2 ) to be 



n 1 ,n 2 ;t) 



(38) 



p[ e (p,n 1 ,n 2 ;p f ,n f 1 ,n' 2 ) 



[-iWppt - (7 + 7e)/2] pL(p, ni, n 2 ; p', n' 1; ra' 2 ' 



* i(Qi— ui e i)t VI 



[pL(P + ^ki, ni - 1, n 2 ; p', n[,n 2 



-p I 11 (p,n 1 ,n 2 ;p' - Kk^n^ + l,n 2 



-Pi 1 (p,ni ) n 2 ;p' - hk p ,n[,n 2 



+i X * 2 e i ^ + ^ t p I 12 (p,n 1 ,n 2 ; p' - /jk 2 ,ni,n^ + 1), 



(39) 

where the £ argument has been dropped. By carrying out the trace in Eq. (p7|), one obtains 
terms such as, 

Pn(p; p' - ftki) = 5^ ^11 ( p> ni ' ™ 2 ' p ' _ ^ kl ' 711 + 1; ™ 2 )' 

ni,n 2 

so that the equation of motion for p' le (p;p') has a form which is identical to that which 
would have been obtained using classical pump fields, 



Ple(P; P') = - [KM? + Vie + Upp>) + (7 + 7e)/2] p le (p] p) 

-«Xp [pL(p + ^k P ; p) - Pn(p; p - % 



*p — iSpit V I 



[pL(p + ^ k i; p') - pu(p; p' - hk i 



+tX* 2 e- lS ^ t p I l2 (p;p'-hk 2 ). 



(40) 



In terms of the perturbation series solution 

pIM p) = pff( P ; p') + pUp; p'K v * + Pn(p; pV^ 1 '; 

p{ 2 (p; p') = \pf 2 \p; p') + p+ (p; p')e 1 ^ + P r 2 (p; p'K^' 



JSl2t. 



(41a) 
(41b) 
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where pfl (p; p') are independent of Xp and pf 7 '(p;p') are linear in x P , the steady state 



solution for large detuning is 

p'ie(p; p') ~ -r (xaPiaCp; p' - ^2) + x*ph(p; p' - £ki) + Xppffte; p' - Kk p ) 

By making a change of variables, p — p' — > p — p' — Kk p in Eq. fl36|), one gets, 

g— i(kj,-R— Cl p t) 



(42) 



P'le(R^) 




(27r/;) , A , , rfWe i(p - p > R ^(p-p') 

x {%2 (pi(p; p' + ^ k P2) + p| 2 (p-^ k Pi; p + ^12)) + 
X* (pu(p; p + ^ k Pi) + ph(p-^V; p')) + xpp??(p; p')} 



(43) 



One must now obtain equations for pd(p; p') and p[ 2 {p', p') using the effective Hamiltonian 
Htot= H + V/, solve these equations, and then extract Pi~ 2 (p; p') and Pn(p; p') from these 
solutions using Eqs. (|41a| - |41b| ). 

The terms appearing in Eq. (fHf) may be expressed in terms of the dressed state density 
matrix elements. If one is interested only in terms linear in Xpi one can expand the dressed 
state density matrix to first order in Xp as 



Pip 
Pjj' 



iujiot . 



(44a) 
(44b) 



where I,J— {A , B }, V = {A u Bx}, and J' = {A 2 , B 2 }. The p+(p; p') needed in Eq.(g| 



can be expressed in terms of dressed state density matrix elements using Eqs. (piiP2j[49D as 
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Pii(p;p + fikpi) 



sin 6 l (p)sin6'i(p)p 



sin 6*o (p) cos 6*i(p)p 



cos 9q (p) sin 6*i (p) p 



'b a 1 



+ cos 6*o(p) cos6 l i(p)p 



(45a) 



Pn(p-ftk p i;p) = sin0 2 (p)sin0 Q (p)p^ 2Ao - sin# 2 (p) cos (p)p\ 2Bo 

-cos 2 (p) sin^o(p)PB 2 A + cos #2(p) cos^ (p)Pb 2 b ; ( 45b ) 
P^ 2 (P;P + ^ k p2) = e^ d (-sin^o(p) cos#i(p)p^ ()j4i -sin# (p) sin Oi(v)p\ oBi 

+ cos0 o (p)cos0i(p)pJ oAl + cos0 o (p)sin0 1 (p)pJ oBi ) ; (45c) 
Pi2(p-^ k pi;P + ^12) = e^ d (-sin 2 (p) cos# (p)pji 2j4o -sin0 2 (p) sin 0o(p)Pa 2 b o 

+ cos0 2 (p)cos0 o (p)Pb 2 a ) +cos0 2 (p)sin0 o (p)PB 2 B o ) • (45d) 

The state vector in the Schrodinger representation may be expanded in terms of the 
dressed states for the 0, 1, and 2 manifolds, 



* >= 4 (P)IA) > +c Bo (p)|A, > +c^(p)|A 1 > +c Bl (p)\B 1 > +<f Aa (p)\A 2 > +c s B2 (p)\B 2 > . 



In the following, the momentum labels are suppressed. The Schrodinger equation for the 
dressed state amplitudes is then given by 



and matrix elements of V/^ are given in Eqs. fl33a|433hj) . 

Using Eq. fl46|) along with the ground state decay rate 7 and incoherent pumping to the 
manifold, one finds that density matrix elements for the six dressed states in the three 
manifolds (0), (1), (2), evolve as 



ihc = (H + V Id ) c 



(46) 



where 



H G = diag (E^ , E$ , ,E B \ E^ ,E^ B ) 



(47) 




(48) 
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where 



/ 



\ 



9 A x A x PAxBx PAxAo PAxB a PA ± A 2 PA Y B 2 



PBxAx PBxBx PBxAo PBxBo Pb x a 2 Pb x b 2 



Pd(p,p) = c(p)c t (p') = 



Pa Ax PAqBx Pa a Pa b o Pa a 2 Pa b 2 



(49) 



Pb Ax PbqBx PbqAq Pb b o Pb a 2 Pb b 2 



Pa 2 a x Pa 2 Bx Pa 2 a Pa 2 b Pa 2 a 2 Pa 2 b 2 



y PB 2 Ax Pb 2 b x Pb 2 a Pb 2 b Pb 2 a 2 Pb 2 b 2 J 



and the pumping matrix, A^, has the block diagonal form 



A D — diag(0, A d (p, p'), 0) 



(50) 



where A^(p) has the basic structure given in Eq. modified to allow for p, p' coher- 

ence. In particular, the off diagonal elements of Ad(p, p') that give rise to the nonsecular 
contribution to the line shape are of the form 

Aab(p, p'+Kk 12 ) =A 2 (p+ftki2, p'+Kk 12 ) cos [0 o (p)] sin [0 o (p')]-Ai(p, p') sin [0 o (p)] cos [0 o (p')] • 

The general form of the solution is linked to the incoherent pumping of levels 1 and 
2. For a subrecoil cooled vapor, the pumping rate density for bare state density matrix 
elements Py(p; p') is assumed to be 



where AiV~ l has the dimensions of (volume x time) -1 and can be interpreted as the pumping 
rate to state 1 or 2 in position space. With this form of pumping, p^ (p;p'), must be 
proportional to 5(p — p' — Kk') where k' is some algebraic combination of the pump and 
probe field propagation vectors. To obtain p' le from Eq. fl43|) , one must keep only those 
terms in the integrand of Eq. ( f43[) proportional to 5 (p — p') . The (p; p) term in Eq. 
flUf) makes no contribution to the absorption since it is real and will be ignored from this 



A^p, p') = A^- 1 (27r^) 3 <5(p)5(p - p% 



(51) 



point on. 
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The incoherent pumping of states 1 and 2 populate two different manifolds. The pumping 
of state 1 populates the (p = 0, rii, n 2 ) manifold while the pumping of state 2 populates the 
(p = —Kki2, rii, no) manifold since this manifold involves state 2 with zero momentum. Thus 
in viewing absorption or emission, two distinct initial state manifolds must be included, 
leading to the possibility of sixteen rather than eight components of the spectrum. Here we 
set A 2 = 0. 

Substituting Eqs. (|45aj - |45clD into Eq. (f43"D and using the steady state solutions for pf r 
and pjj, obtained from Eqs. (|44a| - [44bD and Eq. (48|), one finds, after some manipulation, 
the final expression for the absorption coefficient, 



(-) 

V a o/ 

(-) 

\O0jru, 



a \ 



\G\ A 



AB 



a \ 

a 0/ns 

sin 2#o 



ij V sin 6 1 -l cos 9i ) (T AoAl (A') + T BoAl (A')) 



rWA') 
rWA) 

IWAO 
Ta a 2 (A') 

iWA) 

^A B 2 (A') 



^B B 2 (A') 



|A| 7 2 ((4°^ + 7 ^ 
- Ur) cos 0! + l - sin 0^ (r AoBl (A') + T BoBl (A')) 

+ _ v -2 _ 2i , cot 2 # ) ( sin 2 ^ (r A)A2 (A') + T BoA2 (A')) + cos 2 # 2 (r AoS2 (A') + r Bofla (A'))) 

^a Ai(AO; 



A' - U hKl - -(^AB - ^AI?) ) U AB + 7 2 



(°) 

A / 1 / (0) , (1) \ \ (0) 2 



A' . , *- I. ,(0) , X}-) \ \ , W) , .,2 

A - Wfik p i + T^AB _ 10 AB) J AB + 7 



,(!) 



,(o) 



1 

<„\ + 2 ■ 
I 

vi - 2 ' 

1 



A' 4-/, V/ > (0) 4-/ > {2 h ^ , ) (0) 4-o/ 2 

A + Wakp! - ^AB + ^abJ ^ab + 7 



-(A' + ^ + ^ + ^I^b+T 2 



1 



A' + u fatvl - ^(u% - u} y AB ) ) uj k a ' b + 7" 



,( 2 ) 



,(0) 



-^BoAi(A'); 
^AoBi(A'); 
L BoBl (A'); 

^a a 2 (A ); 
-^BoA 2 (A'); 

-^AoB 2 (A'); 
Lb 0j b 2 (A'); 



where Aab = srn (2#o) an d ( — ) is given by Eq. fl27a| ). The reason the absorption 

\ a ° / sec 

coefficient is expressible as a sum of the secular term plus a non-secular term is linked 
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to the fact that the secular approximation consists solely of neglecting the off-diagonal 
components of A<j. Since the first order solutions, p\ v and Pjji, are linear in the pumping 
terms, ( ^- ) contain terms proportional to and A# while [^-) is proportional to 

\ a ° / sec V a ° / ns 

Aab- This simplification would not occur for a more complex decay scheme for states |1 > 
and 1 2 > since the decay would couple density matrix elements in a field dependent manner 
(see [H). 

The nonsecular term, (j^ j , consist of dispersion-like structures centered at the same 
locations as the resonances in ( — ) .In the secular limit, ( — ) < f— ) and ( — ) 

Wsec ' \ a °Jns \ a °Jsec \ a ° J ns 

can usually be ignored. Notice that if one chooses ip and rj such that a pair of absorption 
lines in (j^j vanish, then the corresponding terms in (^^j a l so vanish so that (^^j is 
identically zero. However, this will not be true for the gain terms in ( ^- ) since the cor- 

V a ° / sec 

responding terms in (^^j have a different interference coefficient, {rf — r^ 2 — 2ip cot 26$). 
Figure 6 shows a plot of the non-secular absorption coefficient for the same parameters as 
Fig. 3. In this plot, the non-secular terms are ~ 1000 times smaller than the secular terms. 
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Figure 1. Schematic diagram of atom-field system. 

Figure 2. (a) Transitions between the manifold and the 1 and 2 manifolds leading to 
probe gain or absorption in the bare state picture. The states in the 1 and 2 manifolds 
are displaced from the states in the manifold by an amount ~ hQ±. The Rabi frequencies 
shown are those which couple the states from the manifold to the 1 and 2 manifolds, (b) 
Energy levels in the dressed state basis. The center of the 1 manifold has an energy huio 
above the center of the manifold and similarly, the center of the 2 manifold is hu 2 o below 
the manifold, (c) Illustration of the coupling of the pump and probe fields to the dressed 
states for the \A >^ \A\ > transition. The Rabi frequencies shown are those that couple 
the bare state |1 > component of \A > to the bare state components of \Ax >. 

Figure 3. Plot of ( f- ) for A 2 = 0, A1/7 = 1, ip = -1 and 77 = 2 showing all eight 

V a ° / sec 

absorption and emission lines. The detuning is #12/7 = 300 and the two-photon pump 
Rabi frequency is |G|/7 = 250. The recoil energies are uj^ri/l = 40, u)Rk. pl h = 60, and 
Kk pl ■ k 12 /M7 = 80. Note that 7 7 e where 7 e is the excited state decay rate. 

Figure 4. Plot of (j^j showing destructive interference. Parameters are the same as 
Fig. 4 except for ip and rj. The solid line corresponds to ip — — 1 and rj = Vtan^i = 0.8383 
while the dotted line corresponds to ip — +1 and rj = \J cot Q\ = 1.1928. 

Figure 5. Plot showing sixteen lines in the ( — J probe spectrum for A2/7 = A1/7 = 1, 

V a ° / sec 

tp — +1 and 77 = 0.1. The detuning is #12/7 = 500 and the two-photon pump Rabi frequency 
is |G|/7 = 750. The recoil energies are c<; ftkl2 /7 = 50, w?uc pl /7 = 75, and ftk pl -k 12 /M7 = —50. 

Figure 6. Plot of (^^j for the same parameters as Fig. (3). The non-secular absorption 
coefficient has dispersionlike structures at the same location is the line centers of (j^ j 
The amplitudes of these non-secular terms is typically ~ 1000 times smaller than the secular 
line strengths, consistent with ^/oj^b = 0.00177. 
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